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I. INTRODUCTION 

In all early Hartree-Fock-Slater (HFS) calculations^’^ of radiation- 
less transition rates in atoms with an inner-shell hole, frozen orbitals 
have been assumed. The wave function of the continuum electron is 
calculated in the V potential of the initial configuration, as though 
the continuum electron sees the same potential as the bound electrons 
in the initial state. In this approximation, tfie exchange potential used 
in finding the continuum wave function is the average exchange potential 
of the bound electrons, and not the exchange potential between bound and 
continuum electrons. 

Calculations chat include the frozen-orbital approximation have been 

reasonably successful in predicting probabilities of energetic Auger 
1 2 

transitions. ’ For Coster-Kronig and super Coster-Kronig transitions, 
on the other hand, in which the continuum-electron energy is usually quite 
low, large discrepancies are found between existing theoretical results 
and measured rates. It is of interest to examine to what extent these 
perplexing discrepancies can be ascribed to the neglect of exchange and 
relaxation effects. Exchange has indeed been found to be significant in 

p 

low-energy electron scattering from atoms and in photoionization near 

Q 

threshold, and its influence on low-energy radiationless transitions 
has recently been noted by McGuire,^ In this paper, we report on an 
investigation of exchange and relaxation effects in low-energy Coster- 
Kronig and super Coster-Kronig transitions, applied specifically to 
transitions that deexcite the [3p] vacancy state of atomic Zn. 
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II. THEORY 
2-fi 

In frozen-orbital HFS calculations of the Auger transition 
probability, the wave function of the continuum electron Is calculated 
In the same potential as that of the bound electrons. The radial equation 
then Is 

[■ i * ^ * '"’■*] 

where the potential Is 

V(r) = 2 q(nji)'^g(n£n)i|r) -6[^^lp(r)|| . (2) 

The quantities {n£n'i,'[r) are deflned^^ as 

■^l^(njn'ii'lr) = P„,(r) ^P„,j,(t)dt +1 P„,(t)-{4r Pn',.' 

'■ (3) 

Here, p{r) Is the spherically averaged total electronic charge density. 

The Slater exchange potential [the last term in V{r)] is the "average" 
exchange potential for the bound electrons and does not represent the 
exchange interaction between bound and continuum elections. The Slater 
exchange term introduces an attractive potential that pulls the continuum 
electron inward. The effect of this term on super Coster-Kronig rates is 
discussed in Sec. V. 

We now proceed to derive the Hartree-Fock (HF) equation for the 
continuum electron which contains the correct exchange interaction between 
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bound and continuum electrons. 

The wave functions of the Initial and final states of the atom are 
assumed to be represented by linear combinations of Slater determinants 
with normalized spin orbitals. In LS coupling, the wave functions are 
characterized by the total orbital (L) and total spin (S) quantum numbers. 
For the specific examples considered here, in which super Coster-Kronig 
and Coster-Kronig processes fill a 3p-shell vacancy in atomic Zn, the 
transitions can be described as follows: 


ls^2s^2pSs^3p®3d^°4s^ V>ls^2s^2p^3s^3pSd°(S(,L(,)ca' ^P, 
and 


2.'=! ,3,5 


ls^2s^2pSs^3p^3d^°4s^ ^P->ls^2s^2p®3s%®3d^4s^S(.'Lj.' )e£' ^P, £'=1,3, 

Here, e£' denotes the energy and orbital angular momentum of the Auger 
electron in the continuum. 

We calculate the initial- and final -state wave functions by the 
HF approach. We assume that the core of the final -state atom is unaffected 
by the e£' continuum electron. For the bound electrons, the HF equations 
for the radial wave functions can be constructed following the standard 

n 12 

procedure. * The HF equation for the continuum state is obtained 

13 

following Seaton's method. 

The final state is formed by adding an electron to the unperturbed 
atom with N bound electrons and nuclear charge Z. The antisymmetric 
(N+1 )-electron wave function can be expanded in terms of basis functions 
which are completely antisymmetric under interchange of the coordinates of 
a pair of bound electrons: 
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N+1 

4'(ijlX) = \ ^ ( -1 (N+1 I X'’ ,x^ ) ; 

'Kr|X \x^) » ^(rlx \r^c^) ' 


(4a) 

(4b) 


.4 

Here, X denotes all the coordinates of the N+1 electrons, except those 
of the electron; Fp(r) is the radial wave function of the continuum 
electron, and r^, o. are the angular and spin coordinates of the i^^ electron, 
respectively. The complete set of quantum numbers required to specify 
the atomic system in the state J system of 

bound atomic electrons is characterized by ctjLjSj, and the continuum 
electron, by 

u J 

13 

The HF equation for the continuum state can be obtained from 

X dX-<”*'>df|,^,do„+, = 0. (5) 

Upon integrating over angular coordinates and summing over spin, we find 
for the radial equation of the continuum electron 


i- 

^ ^^N+l 


"^N+l 


2Z 

"n+1 


+ k‘ 


] 


'"r^'^N+l^ 




- 2Wj,p(r^^^) 


a 

where A., is the off-diagonal Lagrange multiplier, and we have 


( 6 ) 
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Vrr(^N+l) = ‘‘'*(^’1^ 




-(N+l)i 


N+rN+1 


r(N+l); 


N+l,i 


N+1^N+1' 


^rr^'"N+l^ 




X ^^N+l^®N+] * 

The direct term V^j. and the exchange term Wj,j, were evaluated by the 

method of Hartree^^ and Racah,^^ and the results were checked against other 

15 16 

existing calculations. ’ 

2 

We shall work out the specific cases of 3p-3d ct'{£'=l ,3,5) and 
3p-3d4se£' {t'=l,3) radiationless transitions in atomic Zn. For these 
transitions, Vpp and can be expressed as follows. 

''rr ^ a*^2<3^3d|r) + b*^(3d3djr) (7) 


Wj,j, =-[2(2£'+1)]'*^^^ q(ns)*^,(nsEt'lr)P(ns) 
ns ^ 

- y^|[3t'/(2£’-H)(2£'-l)]^^,.,(npc£'|r)P(np) 

np 

+C(3£'+l)/(2£'+3){2Jl'+l)]'^j^,^^(npe£’ |r)P(np)J 
-cl^._2(3de£'|r)P(3d) - d‘^j^,(3dc£'|r)P(3d) 

,(4se£' |r)P(4s). 
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Here, the sum in Vj.^, extends over all occupied shells (i.e., Is, 2s, 

2p, in the final atom, after the Auger transition. The number of 
electrons in a shell ni, is denoted by q(nji). We have ns = Is, 2s, 3s, 

4s and np = 2p, 3p. The radial wave function of the r\i bound electron 
in the final atom is P(n£). The coefficients a, b, c, d, e, f are given 
in Table I. 

In deriving Eqs. (6), (7), and (8), it has been assumed that the 
orbitals are mutually orthogonal. This assumption involves no loss of 
generality for the continuum wave functions, because there is no open 
shell in the final atom having the same orbital symmetry as the continuum 
electron. In Eq. (3), terms with off-diagonal Lagrange multipliers are 
added in the HF equation for the cp continuum electron to insure ortho- 
gonality. 

Using the orthogonality requirement, one can obtain the off-diagonal 
1 7 

Lagrange multipliers a^_ _ and „ from the HF equation for the 

, t,p jp , I p 

continuum electron [Eq. (3)] and the HF equations for the 2p and 3p bound 

electrons of the final atom. The expression for „ is 

np,ep 


^np,ep a'R^(np3d;3dep) + b'R3(np3d;3dep) 

+ c'R2(np3d;ep3d) + d' R^ (n ;4s ;4sep) , 
where we have np = 2p or 3p, and 

^'^1*^1 '^2^2 ’*^3^3 *^4 4^ " jj^n^ ^^1 ^'"n2£2^’^2^~FiT 


(9) 




is the generalized Slater integral. The coefficients a', b', c', and d' are 
listed in Table II. 
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III. THE RADIATIONLESS TRANSITION PROBABILITY 
18 

From Nentzel's ansatz. the radiationless transition probability 
is given by the familiar formula of perturbation theory, 

CO) 

where the continuum wave function has been normalized to represent one 

19 20 

electron ejected per unit time. * To take into account the effect 

of relaxation, we calculate the Initial and final bound states separately 

in accordance with their respective hole configurations. The lack of 

orthogonality between initial and final electron orbitals with the same 

orbital symmetry causes the expression for the Auger transition probability 

21 22 

to contain more terms than in the usual frozen-orbital theory. ’ For 

p 

the 3p-3d e?. and 3p-3d4sfcA transitions in atomic Zn, however, there are 
no contributions from the one-electron operator because of the orthogonality 
of the angular parts of the wave functions. 

The radiationless transition rates for Zn 3p-shell super Coster-Kronig 

p 

(3p-3d e£) and Coster-Kronig (3p-3d4set) transitions, including the effect 
of relaxation, are given by the following expressions. 

W,i(3p-.3d^) = 4T2|K,<^3p'a|J-|3d^ .[<2p'El|^|3d^> < 3p'|2pJ> 

! 2 

X (quantity of order one) + ... , / 1 ,n 
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wlth 


K, = {ls'|ls)^<2s'|i»s)^<2p'l2p>®<3s’|3s>^<3p'j3p>^ <3d'|3d^^ <4s'|4s^^ 
+ smaller terms containing two or more exchange overlap integrals. 

(3p t-3d4s) = <fi’’^|K2^3p'ci(,|™- 3d4s^ - [<2P 'a|^[3d4s> 

12 12 


X <(3p’ |2p^ + {3p 


12 


^3p'cii 


12 


3d4s)(En|3p^6^^ 

3d3s^^3s ' !4s^Jx {quantity of order one) -r 


( 12 ) 


with 

Kg =^ls'!ls)^<2s'|2s)^<2p'j2p^® <3s' j3s>2 <3p' |3p>^<3d' |3d)^<4s' |4s^ 

+ smaller terms containing two or more exchange overlap integrals, 

In the present calculations we neglect terms with one or more exchange 
overlap integrals, such as ^3p' |2p^ , ^c£{3p^ . Then the radiationless 
transition rates become 


w^.(3p->3d^) =/«‘^K^^k3p 



2 


and 


w^.(3p'^3d4s) = <(3 p’eJI ^ 3d4s^ 


The two-electron Auger matrix elements, expressed in terms of 

generalized Slater integrals, were calculated using Racah algebra, 

21 22 

The relevant expressions * are given in Table III. 


(13) 


(14) 


IV. NUMERICAL METHOD 


The initial and final bound-state wave functions were calculated 

1 2 

with Froese-Fischer's Hartree-Fock program. An iterative procedure was 


-lo- 


used to solve the HF equations for the continuum state, the first solution 

being obtained by Ignoring the exchange terms. The iterations were continued 

until the maximum relative error In the wave functions became less than 
-4 

10 . In general, the solutions converged to this extent within 6 

Iterations for ef and eh states, and within 20 Iterations for cp 

states. This convergence criterion was found to lead to consistency 

to four significant figures In the Auger radial matrix element. 

The 3p-shell super Coster-Kronig and Coster-Kroniq transition energies 

1 2 

were calculated as well with the Froese-FIscher HF program, Including 
relaxation. The calculated energies are compared with experimental results 
In Table IV. 


V. RESULTS AND DISCUSSION 

A. Frozen-orbital approximation 

In order to study the sensitivity of the super Coster-Kronig 

4 2T 

transition probability on wave functions and energy, ’ we have calculated 

these transition probabilities on the basis of several frozen-orbital models. 

To distinguish various models, we shall use N to denote the number of 

bound electrons in a neutral atom (I.e., N»Z). Results for the Zn 3p- 
2 3 1 

3d ( F; G)ef transitions are shown In Fig. 1. To calculate the rates 
N-1 

labeled "HF V ", Initial -configuration HF bound-state wave functions were 
used for both Initial and final bound states; the continuum wave function 
was calculated in the V potential of the (singly ionized) Initial 
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N 

configuration, without exchange. The "HF V “ rates were calculated by 

a similar procedure, the only difference being that the HF bound-state 

N-l 

wave functions of the neutral atom were used. The "HFS V " rates were 

24 

calculated with Herman-Skillman HFS wave functions for the initial hole 

state, employed for both the initial and final bound states; the continuum 

wave function here is the solution of Eq. (1) with V(r) [Eq. (2)] taken 

to be the average potential for the initial state. 

It is clear from Fig. 1 that tne 7n 3p-3d^{^F;^G)ef transitions 

are quite sensitive to the transition energy and to the details of the 

N-1 

wave functions. The frozen HFS V rates, for example, are much larger 
than those derived for the HF model. Most of this difference can be 
attributed to the fact that V(r) [Eq. {!)] in the HFS model contains 
the wrong exchange term. 

B. Exchange and relaxation effects on 3p-shell Coster-Kronig and super 
Coster-Kronig transitions in Zn 

To investigate how exchange affects the transition rates, calculations 

were performed with different continuum wave functions. The continuum wave 

functions were found by solving Eq. {6} (D neglecting the exchange term 

and (2) including the exchange term W^j,. In Figs. 2 and 3 we have 

plotted the continuum wave functions for 3p-3d ( G)ef and 3p-3d ( S)ep 

transitions, respectively; the transition energy was taken to be 42 eV 

25 

in both cases, equal to the average observed energy. Also shown in Figs. 

2 and 3 is the HF 3d wave function of the final atom. 
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It Can be noted In Fig. 2 that there is very little cancellation 
in the overlap between 3d and cf wave functions. The 3p-3d^(^F)ff and 
3p-3d ( G)t;f rates are therefore large and very sensitive to energy and 
wave functions. By contrast, Fig. 3 shows that in 3p-3d^(^S;^P;^D)up 
transitions strong cancellation occurs In the overlap between 3d and rp 
wave functions. These transition rates are therefore small. The inclusion 
of exchange is seen to pull the continuum wave function In toward the origin, 
causing the 3d cf transition probabilities to be increased by 6% (Fig. 2), 
while the 3d^(^S)cp rate is reduced by - 1 % (Fig. 3). The effect of 
exchange becomes less pronounced as the Auger energy increases, 

Using HF instead of HFS bound-state wave functions in the HF continu'.rr; 
equation and in Auger transition probability calculations redujes the 
transition rates by only -3% fore= 42 eV (Fig. 4). 

Relaxation is found to affect the transition rates substantially, 
Relaxation is taken into account by discarding the frozen-orbital approxi- 
mation and using HF bound-state wave functions corresponding to the appropriate 
(different) initial and final states in the calculations. The exchange 
overlap integrals, such as (2p'j3p^ , ^ep|3p^ , In Eqs. (11) and (12) are 
all found to be s 10“^, so that terms involving these integrals could 
justifiably be neglected. For the overlap correction factors we find 

K^*0.9239 and In Fig. 4, the 3p-3d^(^G)ef transition probabilities 

N-2 

are shown, as calculated with continuum wave functions in the V potential 
and with different bound-state wave functions. Relaxation is seen to reduce 
these transition rates by -20%. 
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Also Included in Fig, 4 life results (labeled "HFS that were 

obtained with HFS Initial -configuration bound-state wave functions and 
with a continuum wave function that Is the solution of Eq. (1) with V(r) 
[Eq. (2)] being the average potential of the final two-hole configuration. 
Once again, these HFS results are much larger than the rates calculated 
by other nuthods, due to the error In the exchange term In V(r) of Eg. (2). 

C. Zn width 

The 3p-level widths of Zn, calculated with different wave functions, 
Ci'e listed In Table V. Average experimental Auger energies of Zn vapor^^ 
were used in the calculations. The experimental Mg g level wldth^^ of 
atomic Zn Is listed for comparison. The agreement between the HP 
frozen-orbital result with experiment is probably accidental, Results 
from the relaxed HF model, which can be considered more realistic, agree 
better with experiment than HFS results but still are 50/5 too large. 

Transition rates to various final states, calculated with the 
relaxed HF model, are Indicated In Table VI, 

VI. CONCLUSION 

In HFS calculations of radiationless transition rates, the assumption 
that the continuum electron sees the same potential as the bound electrons 
can cause a large error in low-energy Coster-Kronig and super Coster-Kronig 
cases. The HF approach including relaxation removes some but not all of 
the discrepancy with experiment. 
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The effect of configuration Interaction among discrete states on 
3p-shell Coster-Kronig and super Coster-Kronig transitions of Zn 

4 

has been found to be small. The Importance of the final -state channel 
interactions has recently been pointed out by Howat et al. Calculations 
to include the channel interactions are In progress. 
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TA8LE I. Coefficients ...» f In Eqs. (7) and (8). 


Final hole state 


Configuration 

Term 

i' 

a 

b 

c 

d 

e 

f 

3d^ 


P 

0 

0 

0 

8/15 

12/35 

0 

3d^ 

3p 

P 

-1/5 

0 

0 

-7/30 

6/35 

0 

3d^ 


P 

3/35 

0 

0 

13/30 

102/245 

0 

3d^ 


f 

24/245 

-8/147 

2/5 

27/451 

50/231 

0 

3d^ 


f 

-3/35 

1/21 

-6/35 

2/105 

50/231 

0 

3d^ 

^6 

f 

-55/147 

-13/147 

-12/35 

46/245 

50/231 

0 

3d^ 


h 

-8/21 

-2/21 

-47/380 

50/429 

21/143 

0 

’3d4s 


P 

-1/5 

0 

0 

1/3 

3/7 

0 

3d4s 


P 

-1/5 

0 

0 

-1/3 

3/7 

-1/3 
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TABLE II, Coefficients a', d‘ in Eq. (9). 



Final hole state 
Configuration Term 

a' 

b' 

c* 

d' 


3d^ 

's 

0 

0 

0 

0 


3p 

-23/15 

-12/35 

2/5 

0 


'd 

-1/5 

36/245 

-6/35 

0 

3d4s 


-8/15 

3/35 

2/5 

0 



-28/15 

3/35 

2/5 

-2/3 
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Wsefc' radiationless transition 


lerms of radial Integrals. 


Transition probability 

(3pcp3d3d) + |- (3pEp3d3d)P 
R^ (3pcp3d3d) + j R3(3ptp3d3d)|^ 

j R^(3pep3d3d) + |-R3 (3pep3d3d)|^ 

^1 ^ R^(3pef3d3d) + ^ R 3 (3pef3d3d)i^ 
^1- ~ R^(3pef3d3d) + ^ R 3 ( 3 pef 3 d 3 d) 

3 R^(3pcf3d3d) + | R3(3pcf3d3d) 

5 |^|R 3 ( 3 peh 3 d 3 d )|2 

^|R^(3pep3d4s) + R^(3pep4s3d)|^ 
j|R^ (3pep3d4s) - R^(3pEp4s3d) ^ 

|||- R3(3pef3d4s) + j R^ (3pEf4s3d) 

|l- j- R3(3pef3d4s) + ]■ R^ (3pEf4s3d) | ^ 
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TABLE IV. The 3p-shell super Coster-Kronig and Coster-Kronig 
transition energies of atomic Zn. 
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TABLE V. Zn Mg ^3 level width 


Method Width (eV) 


Theory: HFS 

3.25 

Theory; HFS 

4.49 

Theory: HF 

1.83 

Theory: HF 

1.20 

Theory: Relaxed HF 

2.93 

Experiment® 

2 . 1 + 0.2 


®Ref. 25. 
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TABLE VI. Theoretical Zn 3p-shell super Coster-Kronig and Coster-Kronig 
transition probabilities according to the relaxed HP nx)del. 


Transition 

Final term 

Transition rate 
{mini atomic units) 

3p-3d3dep 


0.838 



0.575 



0.740 

3p-3d3def 


8.389 



26.637 


^6 

61.370 

3p-3d3dch 


0.021 

3p-3d4scp 


4.690 


3d 

3.651 

3p-3d4scf 


0.606 


0.048 
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Figure Captions 

FIG. 1. Probabilities of the 3p-3d^(^F;^6)Ef transitions In atomic 
Zn, according to frozen-orbltal calculations with various wave functions and 
potentials; Hartree-Fock in the singly Ionized (V ) and neutral -atom 
(V^) potential, and Hartree-Fock-Slater in the potential. To illustrate 
energy dependence of these rates, the transition energy is treated as a 
variable. The measured average transition energy is 42 eV (Ref. 25). 

FIG. 2. The Hartree-Fock bound-state 3d radial wave function of atomic 
Zn and the cf HF continuum wave function, calculated in the potential of the 
final atomic state, with and without the exchange term in Eq. (6). 

FIG. 3. The HF 3d bound-state radial wave function of Zn and the tp 
HF continuum wave function, calculated with and without the exchange term 
W^j, in Eq. (6). 

2 1 

FIG, 4. Theoretical 3p-3d ( G)ef radiationless transition probabilities, 

as a function of Auger energy, calculated according to various models. Relaxed 

Hartree-Fock, frozen Hartree-Fock, and frozen Hartree-Fock-Slater bound-state 

wave functions are used, all three with a Hartree-Fock continuum wave function 
N-? 

in the V potential of the doubly ionized final-state atom, with and without 
exchange between continuum and bound electrons (Eq. 6). Only for the top 
curve (diamonds), the continuum electron is described by a Hartree-Fock-Slater 
wave function [Eqs. (1) and (2)]. 
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